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Abstract 

In this paper, we introduce a numerical renormalization group procedure 
which permits long-time simulations of vortex dynamics and coalescence in a 
2D turbulent decaying fluid. The number of vortices decreases as N ~ i 
with £ ~ 1 instead of the value £ = 4/3 predicted by a naive kinetic theory. 
For short time, we find an effective exponent £ ~ 0.7 consistent with previous 
simulations and experiments. We show that the mean square displacement 
of surviving vortices grows as (x 2 ) ~ t 1+ ^ 2 . Introducing effective dynamics 
for two-body and three-body collisions, we justify that only the latter become 
relevant at small vortex area coverage. A kinetic theory consistent with this 
mechanism leads to £ = 1. We find that the theoretical relations between 
kinetic parameters are all in good agreement with experiments. 



PACS numbers: 47.10.+g, 47.27.-i 
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I. INTRODUCTION 



In recent years, a great deal of work has been devoted to the study of two-dimensional 
turbulence. Two-dimensional turbulence is not only relevant to the study of geophysical and 
astrophysical flows, but it is also far more accessible to modern computers and experiments, 
since the measurement and the visualization of the velocity and vorticity fields are much 
easier than in D = 3. In addition, two-dimensional turbulence has deep connections with 
other fields of physics such as electron plasmas in magnetic field PJ and stellar dynamics [0 . 

For the specific problem of two-dimensional decaying turbulence, recent experimental 
|3|-[5|] and theoretical p|-|TE[ works have emphasized the importance of coherent vortex dy- 
namics during the fluid decay. This process essentially consists in three stages: during an 
initial transient period, the fluid self-organizes and a network of coherent vortices appears. 
Once the coherent vortices have emerged, vortices disappear through mergings of like-sign 
vortices, such that their number N decreases and their average radius r increases, in a 
process somewhat reminiscent of a coarsening dynamics |16[]. During this process, and in 
the limit of small viscosity, energy remains constant. When only one dipole (or very few) 
remains, it finally decays diffusively. 

From the theoretical point of view, the "coarsening" stage is certainly the most interesting 
as, in principle, it can extend on an arbitrary long time period. In this regime, the main 
question arising concerns the existence of universal features including the decay exponent £ 
(N ~ i ), and other exponents which describe the time evolution of quantities such as the 
average vortex radius, the enstrophy or the kurtosis. 

In this paper, we first describe an effective model for the vortex dynamics and review the 
main experimental and numerical results concerning the temporal evolution of the physical 
quantities listed above. In section II, we show that the surviving vortices have a hyperdif- 
fusive motion with an effective diffusion coefficient D ~ t^ 2 , and a flight time distribution 
behaving as -P(t) ~ r^ 2-3 . We then consider a "naive" kinetic theory for the vortex decay 
dynamics predicting £ = 4/3. In section III, we introduce a numerical renormalization 
group (RG) procedure which permits very long simulation times. Although the numerical 
diffusion coefficient is well described by the preceding kinetic theory, the decay exponent is 
found to be significantly lower than expected (£ ~ 1 instead of £ = 4/3). In section IV, we 
derive an effective dynamics for 2 and for 3 neighboring vortices subjected to the effective 
noise due to far away vortices. Within these simple models, we relate the average merging 
time to the decay exponent £ found in the RG simulations. Our main conclusion is that 
the lower than expected value for £ could be explained by the fact that two-body collisions 
are irrelevant at large time, whereas three-body collisions predominate. In section V, we 
present a simple kinetic theory taking these three-body collisions into account and yielding 
£ = 1. The importance of three-body collisions in vortex dynamics was previously pointed 
out by Novikov flT7| ], in a different context. Throughout this paper, we compare our results 



to recent experiments M| and find a very good overall agreement. 



II. KIRCHHOFF MODEL 
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A. Generalities 



As we are mainly concerned with the coherent vortex dynamics and merging processes, 
it is natural to focus on the effective behavior of the sole vortices, neglecting the incoherent 
background. 



The route to such an effective model starts with the work of Kirchhoff |L8j who obtained 
the equations of motion of point-like vortices in the zero viscosity limit. Vortices follow a 
Hamiltonian dynamics where the vortex center coordinates X, and are conjugate variables: 



r ; 



dxi 
~dt 

r 

*dt 



dH 

dH 

dxi ' 



where Tj is the circulation of vortex i, and where H denotes the following Hamiltonian: 



# = -£ r,r, in(i 



(i) 

(2) 
(3) 



where rij stands for the distance between vortices i and j. These equations of motion can 
be more explicitly written: 



dX{ ^ s—-\ Uij 

~dt ~ ~ ^ J r^' 



7 K) 



(4) 
(5) 



These equations are strictly valid for point-like vortices and cannot describe vortex mergings. 
This should be accounted by hand by defining ad hoc merging rules as introduced in [|7].|9].|T3| • 
The authors of determined a criterion for the merging of two like-sign vortices of radii 
ri and r 2 , within the elliptical-moment model [HJ. They found that, for r\ < r 2 , collapse 
can be observed for an initial vortex separation d < r c = ar 2 + 6r^/r 2 , where a and b are 
numerical constant of order 1 (a ~ 2.59, b ~ 0.61 ||). Using r\ < r 2 and a > 3b, one easily 
obtains 



a + b 



(ri + r 2 ) <r c < ah + r 2 ) 



(6) 



which shows that r c is of the order of the mean radius. In |13|], the authors in fact used 
r c = a'{ri + r 2 ) (a' ~ 1.7). Thus, if one considers a collection of vortices of the same typical 
size, it is clear that choosing the same critical distance for all mergings, equal to the average 
radius r of the population of vortices, cannot drastically affect the model properties. This 
was actually verified in |16| . Now that the merging criterion has been given, the properties of 
the vortex resulting from the merging of two like-sign vortices must be specified. Motivated 
by experiments and numerical simulations ,jl2 , 13|| , the authors of |12] , ^3|| and 

assumed that the average peak vorticity u is conserved throughout the merging process, 
as well as the energy since the inviscid limit is considered. As the total energy scales as 
E ~ Nu 2 r A (with a possible In A correction), this shows that Ar 4 should be conserved, 
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or equivalently that a vortex of radius r' = (rf + r 4 ) 1 / 4 results from the merging of two 
vortices of radii r± and r 2 . Note that this conservation law is consistent with the observed 
slow enstrophy dissipation [|l^ . 



Numerical simulations of this model starting from a population of vortices having the 
same typical radius results in a narrow radius distribution at all subsequent times pl,|T6 . 
Moreover, it is observed that the number of vortices decays as a power law 

N(t) (7) 

(i + tAo) ? 

with £ pa 0.70 — 0.75, much smaller than the exponent predicted by Batchelor theory (£ = 2) 
p0|| . The conservation of the total energy and mean peak vorticity leads to the occurrence 
of only one independent exponent for the time evolution of physical quantities 



N(t) ~ R- 2 ~ r*, r ~ (8) 
Z ~ t~V 2 , K ~ t«/ 2 , (9) 

where i? is the typical distance between vortices, and Z and if are respectively the enstrophy 
and the kurtosis. 

The exponent £ and the predicted scaling laws are consistent with experiments |3|-[5| 
and direct numerical simulations || of Navier-Stokes (NS) equation (using a hyperviscous 
dissipation term). 



B. Limitation of numerical simulations 

The Kirchhoff simulations and actual experiments cited above were only carried out for 
very short time. In |9j (as well as in experiments the number of vortices decays 

by less than a factor 4 at the maximal accessible time t max - Raw data show significant 
curvature on a log-log plot, hence the introduction of an extra fitting parameter to i n 
(see Eq. (|7|)). Since to ~ ^max/3, the simulation time is of the order of the transient time to, 
and the scaling regime for t ^> to is probably not reached. As a matter of fact, the actual 
exponent £ obtained by measuring the logarithmic slope at the final time is of order £ pa 0.6, 
as obtained by Benzi and co-workers ( see 8 ^ so section III.B) 



It would be interesting to explore the domain of lower vortex density n, as in all sim- 
ulations performed so far the mean free path (of the order of the typical distance between 
vortices [11,^1]) remains of the same order as the radius size. In other words, the fraction 
of occupied by the vortices, 

Nnr 2 2 fr\ 2 . , 

s = -^ = nvrr = ^ l^-) , (10) 

remains quite large in the early time of the dynamics. Because of the scaling laws of Eq. (§), 
R grows faster than the mean radius r, such that vortices become effectively more and more 
point-like, a regime which seems to be out of numerical reach, and which should develop for 
t > t . 
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C. Kinetic theory for the Kirchhoff model 



1. Diffusion coefficient 

In the absence of mergings, the chaotic Kirchhoff dynamics is known to lead to an effective 
diffusive motion of the vortices |II^ , [lT| , |2~T| | . The diffusion coefficient can be calculated by 
computing the fluctuation time T(v) for a given vortex velocity v and averaging the quantity 
v 2 T(v)/A over the velocity distribution. This calculation has been extensively described in 
plfl and we present here a simple heuristic argument leading more directly to the same 
result. 

Using Eq. ([|) and Eq. (^), the average velocity squared is 

3 y 

where we have neglected the contribution of off-diagonal terms obtained when squaring 
Eq. (£|) and Eq. @. If we assume the vortices to be uniformly distributed on average, and 
the circulations to be equal up to their sign, we then obtain 

„ NT 2 ^- f L ^ ~ 2nT 2 \n(L/r) ~ nT 2 In N, (12) 

Li Jr X 

where we have introduced the vortex typical radius r as a natural cut-off. This expression 
already obtained in has been qualitatively checked in [IT| and is confirmed by our 



simulations of section III.B. 

It is then natural to assume that the mean free path / is of order R, the typical distance 



between vortices, as proved in [pijj . We then get the expressions for the mean free time r 
and the diffusion coefficient: 

/ I 2 
I ~ R, r (nTVh^N)- 1 , D rVlniV, (13) 

V T 



in agreement with our more sophisticated treatment [[21 



Now, if we include merging events, these different quantities are expected to vary with 
time as both the density n and the typical circulation T do. If we drop logarithmic corrections 
for now, we obtain 

I rs, t €/2 , r ~ t f/2 , v ~ const., D ~ t f/2 . (14) 

Note that this expression for the diffusion coefficient D differs from that obtained in || 
(D ~ t 3 ^ 4 ). Indeed, the authors of 0] used the merging time r mergi to compute D, instead 



of the fluctuation or mean free time which is relevant here [|2T]JTT|1 . A naive estimate of this 
merging time is addressed in the next subsection. 

Finally, we predict that the mean square displacement of surviving vortices (or test 
particles) in the decaying fluid should behave as pll 



(x 2 )~r, with i/ = l + |. (15) 
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In U, the authors found £ ~ 0.7, which leads to v ~ 1.35, using Eq. ([15]). This must be 
compared with the experimental value v ~ 1.3 for vortices and i/ w 1.4 for test particles 
moving along the current lines of the fluid. Eq. (Tf5| ) is also in good agreement with our 
simulations of section III.B. 

Note that this hyper diffusive behavior can be interpreted by invoking a power-law de- 
creasing flight time distribution (time between two deflections or large velocity fluctuations). 
Let us assume that this distribution behaves as 

P(t) ~ r-f, for r -> +oo, (16) 

with fi > 2, such that the average fluctuation time exists. Then, after a time t = Y^L\ T % ~ 
m(r) (m deflections), and using that (v 2 ) is essentially constant, the mean square displace- 
ment reads 

(x 2 ) ~ (t; 2 ) r 2 \ ~ m J* t 2 ^ dr ~ t 4 "^. (17) 
This shows that within this interpretation v = 4 — //, or 

M=3-|. (18) 

Using ^ w 0.7, we find \i ~ 2.65, in good agreement with the experimental value \i = 2.6±0.2 
measured in ||. 



2. Naive kinetic theory 

The merging time r merg . is the typical time between two merging events involving the 
same vortex. If we assume that vortex mergings occur whenever two like-sign vortices stand 
at a distance less than r c ~ r, a classical cross-section argument leads to 

r merg . ~ (nvr-y 1 . (19) 

If we assume a scaling regime where n(t) decays as a power law, r merg . must behave linearly 
with time [E^| as 

dn n n , , 

* ~ " e 7 ~ (20) 

Using the scaling equations Eq. (H) and Eq. (|14]) , and the above expression for r merg ., we 
finally obtain 



' mcrg. u 



t^~t3« ; (21) 



neglecting logarithmic corrections in v. The constraint r merg . ~ t then leads to £ = 4/3, well 
above the measured value £ ~ 0.7. Note that our argument is fully consistent with a direct 
simulation performed by Trizac p2| , who found £ ~ 0.7, in a balistic system obeying the 
same conservation laws as in the vortex model, but for which the typical velocity decreases 
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as v ~ t~ 0AJ instead of being constant. Using Eq. (pip , we indeed predict a decay exponent 
£ = |(1 — 0.47) ~ 0.707, in perfect agreement with the observed decay exponent. Note that 
the similarity with the value of the decay exponent observed in vortex dynamics is then 
purely incidental. 

In the experiment ||, the authors obtained T merg . ~ £ ' 6 , which strongly indicates that 
they had not yet reached the scaling regime, which is not too surprising as the experiment 
was performed on less than a time decade. As mentioned in section II. B, this should raise 
some doubts about the validity of the apparent exponent £ ~ 0.7. 



III. NUMERICAL RENORMALIZATION GROUP (RG) 
A. Implementation of the numerical RG 

In this section, we address the problem of performing long-time simulations allowing for 
low vortex density and total vortex area coverage. Direct simulations of Kirchhoff equations 
are doomed to failure since each evaluation of a vortex velocity involves the sum over N 
terms (see Eq. (|]) and Eq. ([5])). To access very low densities, one must therefore start out 
with a very large initial number of vortices, which results in a very slow early dynamics. 

The idea of the new approach that we introduce in this section, is to work with a constant 
number of vortices in a reasonable range (N ~ 20 — 100), and progressively increase the 
domain size L by a procedure detailed below. For simplicity, we work with vortices of 
identical radii at all times and thus of equal circulation up to a sign. This is probably a 
reasonable approximation, as previous short-time simulations and experiments have shown 
that the radius and circulation distributions remain narrow at all times M, 
the universal features of the dynamics are expected to survive. Moreover, keeping such 
constant distributions in time should minimize transient time effects due to the fact that 
the scaling distributions are not yet reached (although this problem is partly taken into 
account by the clever procedure introduced in [[]]). 

We thus consider N vortices of radius r and circulation T = ±umr 2 (N/2 vortices of each 
sign) in a box of initial linear size L, such that the initial density is n = N/L 2 . Periodic 
boundary conditions are considered and Kirchhoff equations of motion are adapted to this 
situation |§||. Compared to Eq. (f|) and Eq. ([|), the velocity induced by a vortex j on a 
vortex i is only significantly different for vortices at a distance of order L, so that the physics 
is not modified particularly for vortices at a typical distance L/\/ r N or less. 

Vortices do obey Kirchhoff dynamics until two like-sign vortices meet, i.e. their distance 
is less than r c = 2r (see the discussion of section II. A). Both vortices are merged, and all 
radii and circulations are updated to 

i 

Y r, V = ±uitr' 2 . (22) 



18, 23]. Therefore, 




The density is updated accordingly: 



A new vortex of the same sign as the vortex which just disappeared is then introduced in 
the box at a random position and has radius r' and circulation V. The number of vortices in 
the box is then restored to its initial value N. All the distances are then scaled by a factor 



N 
N-l 



1/2 



This renormalization procedure ensures that the new density takes the correct value 

, N N-1N N-l 
n= L^ = —V = (25) 

and that the quantity n'r = nr 4 is conserved, ensuring the conservation of the energy per 
area unit. If N is large enough, one may expect that the introduction of a new uncorrelated 
vortex after each merging should not affect the dynamics, especially at large times for 
which the merging time is much larger than the mean free time. In this regime, the newly 
introduced vortex which has only a probability iV -1 of being involved in the next collision 
has plenty of time to get "randomized" as r merg . 3> r. 



B. Numerical results 

We have performed long-time RG simulations with N = 10, 20, 40, 60, 80 vortices, reach- 
ing final densities as low as 2 x 10 _4 n . Except for the case N = 10 which seems to decay 
faster, the different density plots are essentially independent from the actual number of 
vortices involved in the RG (see fig. 1). The long-time decay exponent is estimated to be 
£ = 0.99 ± 0.01 significantly higher than the expected value £ ~ 0.7 but still well below the 
naive estimate £ = 4/3 obtained in section II. C. 

We have also measured the average mean square displacement of surviving vortices. 
The motion is found to be hyperdiffusive with a diffusion exponent v consistent with the 
prediction of section II. C and |2l| . Indeed, we find v = 1.50 ± 0.01, to be compared with 
v = 3/2, admitting the value £ = 1. This is illustrated in fig. 2. Note that according to 
section II. C.2, we thus predict a flight time distribution behaving as P(r) ~ r~ 5 / 2 , for large 
r. 

In the range n/n® = 0.8 — 0.2, corresponding to the range of density obtained in previous 
experiments and simulations |6|,|9|, p!3|Jl6| , we indeed obtain an apparent exponent of 

order £ « 0.7. In fig. 3, we compare the direct simulations of P,|T6| where vortices were 
allowed to develop a scaling radius distribution. After fitting an arbitrary time scale, it 
appears that our RG simulations are in good agreement with these previous works, although 
our RG simulations extend to almost three more decades in time. 

We have also measured the mean square vortex velocity which is expected to behave as 

(v 2 ) ~ nT 2 In 

This behavior is confirmed by our RG numerical simulations as shown in fig. 4. 



m 

At), 



ln(t). 



(26) 



8 



We have also performed RG simulations introducing a distance cut-off in Kirchhoff equa- 
tions Eq. (^) and Eq. @, replacing r?- by (rfj+r 2 ). Indeed, although like-sign vortices cannot 
approach each other closer than a distance 2r (otherwise they merge), opposite-sign vortices 
can, which is quite unphysical as it generates very fast traveling dipoles. Introducing this 
cut-off results in a physical upper cut-off of order T/r for the maximum velocity of these 
dipoles. The number of vortices initially decays slightly more slowly, although the long-time 
decay exponent remains fully compatible with £ « I . 

It is interesting to compare our results to new direct simulations where the radii and 
modulus of circulation are maintained equal for all surviving vortices. Even starting from 
N = 2000, it is hard to reach low densities in a reliable way. As exemplified in fig. 5, these 
direct simulations follow our RG calculations before decaying faster beyond a breaking time 
tjV; as the density approaches the minimum reachable density n/riQ = 2/N. We observe 
that the breakdown occurs sooner for samples with a decreasing initial number of vortices, 
i/v is in fact of the order of the minimum time after which some samples had reached the 
minimum possible density (only 2 opposite vortices left). Still, fig. 5 lends credence to the 
claim that the large N limit finally reproduces the RG results. Note finally that for these 
direct simulations, the actual logarithmic correction in the mean square velocity behaves 
differently from that in the RG simulations as the box size L remains constant. In this case, 
it slowly decreases as 

W-^ln^-ln^), (27) 

up to a time of order t* <~ tjy. Note that this subdominant difference between RG and direct 
simulations could be responsible for a slight discrepancy in the apparent decay exponent 
observed in both kinds of simulations. 



IV. MERGING TIME IN EFFECTIVE TWO-BODY AND THREE-BODY 

DYNAMICS 

A. General form of the merging time 

In this section, we are concerned with an alternative way of evaluating the merging time 
"Tmerg. as a function of the physical parameters n, r, T... On the basis of simple dimensional 
analysis, r merg can be expected to be proportional to the mean free time r multiplied by an 
arbitrary function of the only dimensionless parameter nr 2 . Expecting power-law behaviors 
for these quantities, a natural ansatz is 

T 

Tmerg. ~ ~, (28) 

(nr z ) a 

where a is an exponent to be determined. Dropping logarithmic corrections in r (see 
Eq. (|13D), we obtain 

Tmerg. ~ (nr(nr 2 ) a ) -1 ~ cu" 1 (nr 2 ) . (29) 
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Imposing that r merg . must be proportional to the actual time (see section II. C), and using 
the scaling equations Eq. (|8|), we find the relation between the decay exponent £ and a: 

Note that the naive expression of r merg . obtained in section II. C can also be written in this 
form: 

Tmerg. ~ (rivr)' 1 ~ U^far 2 )' 3 ' 2 , (31) 
corresponding to a = 1/2 and thus £ = 4/3. 



B. Theory of the effective two-body dynamics 

We now consider the effective dynamics of two nearby like-sign vortices, assuming that 
the (N — 2) other vortices are at a distance at least of order R. The velocity induced by the 
other vortices on one of these two vortices can be written as the sum of the velocity induced 
on the center of mass (at r = (r x + r 2 )/2) plus a small correction: 

^ = - E V >% ( 32 ) 

#1,2 ' y 

= v 0x + fa, v + fa, v r,!$C^, (33 ) 

#1,2 %' ^1,2 'Oj 

= w 0x + Sxxrja + Syirjb, (34) 

where we have neglected other corrections involving higher powers of 8x\ = X\ — Xq and 
fyi — Ux~ Ho- After a straightforward calculation, we get a similar equation for v ly : 

v\ y = v 0y + 5x 1 i] b - 5y 1 r] a . (35) 



Note the antisymmetric structure of Eq. ( |3"4D and Eq. fl3"o] ) resulting from the Hamiltonian 
nature of the dynamics. 

We will now assume that Vq x , vq v , r\ a and r\ b can be considered as random Gaussian 
variables. Their second moments are respectively 

<<4> = «> = ^ ~ nY\ (36) 
up to a logarithmic term (see section II. C and | 16| , pT ,2T|), and 

= «> =N f(^~Nf^ f R dx ~ „T». (37) 

We have used i? as the lower cut-off since the other vortices were assumed to be at a distance 
at least of order R. Note that rj a ^ have the dimension of an inverse time and are simply 
of order r _1 ~ nT. It is natural to assume that r] a ^ have the same correlation time as the 
velocity, namely the mean free time r. 
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Assuming that Vo X} y and r] a ^ are Gaussian noises of correlation time r leads us to write 
a simplified effective Langevin equation describing the dynamics of these four quantities: 

du u (u 2 ) , , 

(38) 

where u = Vo xy , r/ a ,b, and w u are independent 5-correlated white noises. Such a Langevin 



equation was recently introduced to describe velocity fluctuations in |pJ]| , with reasonable 
success. 

We are now ready to construct an effective two-body system in order to study the merging 
time. We consider two like-sign vortices of identical circulation in a square box of size R 
with periodic boundary conditions. These vortices are submitted to their mutual advection 
and to the effective noise induced by the other vortices at a distance greater than R. If we 
define x = xl — x2, y = yl — y2 and d = \fx 2 + y 2 , we get from Kirchhoff equations and 
Eq. (HD and Eq. ([§) 

dx y , N 

~dt = d 2 + XVa + VT]bl 
dy x 

~dt = d 2 + XVb ~ y7]ai 

as the average induced velocity cancels out. After expressing time in units of r (t — > t/r) 
and distance in units of R (x — > x/R, y — > y/R), and noting that T ~ i? 2 /r, we end up with 
a dimensionless equation of motion as anticipated in the previous subsection: 

dx y / a t \ 

dt = ~d 2+XVa + yV ^ ( ' 

dy x , 
Tt = j 2 + x V , - y Va , (42) 

where rj a ^ are independent Langevin random variables, of unit mean square average and 
correlation time. 

Both vortices are initially randomly placed in the unit box at a mutual distance d > OAR 
and their relative distance evolves according to Eq. ( |4"1"D and Eq. Q4"^ ) until d becomes smaller 
than the scaled dimensionless parameter d c = 2r/R = 2Vnr 2 , which defines the merging 
time. As anticipated above, the average merging time in units of r can only be a function 
of this parameter d c , leading to Eq. 



C. Absence of strictly two-body collisions 

Numerical simulations of Eq. ( ffID and Eq. ( f42|) lead to the following surprising result: 
both vortices remain at a relative distance greater than a constant G? m j n which slightly de- 
pends on the actual numerical constants in Eq. (|41|) and Eq. (|42|) (to simplify, we have 
assumed the coefficients of (x,y)/d 2 , (rj 2 b ), and the correlation time of rj a ^ to be exactly 
equal to 1). A typical long-time trajectory is shown in fig. 6 and perfectly illustrates the 
absence of collisions when the vortex size is below a certain threshold. Note that if there 
were no noise due to the other vortices, both vortices would strictly remain at the same 
distance, hence producing a circular trajectory. 
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Of course, our result does not prove the absence of collision in the actual iV-body system, 
but strongly suggests that the main assumption according to which all other vortices are at 
a distance greater than R prevents both test vortices from colliding. 

Let us now give a physical interpretation of our result. If both test vortices were at the 
same point, their distance would not vary since the velocity induced by the other vortices 
would be exactly the same for both vortices. Thus, when they are close to each other, the 
effective induced noise is reduced linearly with their distance d as shown by Eq. ( pLTp and 
Eq. fl42|). In addition, as vortices get closer to each other (d ~ r), their relative position 
describes a circle, moving at angular velocity of order 

r 

Q, ~ v r x r _1 ~ — x r~ l ~ uj. (43) 
r 

Because of this fast rotation, the effective fluctuation time of the noise, as seen in the moving 
frame, becomes uo~ l <C r instead of r. Hence, not only the driving noise is reduced due to 
the proximity of the test vortices, but it is also averaged out due to their fast rotation. Such 



a short effective fluctuation time was in fact introduced in Another way of interprating 
the effect of the fast vortex rotation is to note that due to the large difference between the 
system natural frequency uj and that of the noise perturbation (uj' ~ t _1 ), the adiabatic 
theorem ensures that the effective perturbation is reduced by a factor of order exp(— Cujt), 
where C is a constant of order unity ||23|| . 

Conversely, if both vortices start at a distance d < d m i n , this fast rotating pair remains 
stable for a very long time, probably infinite. Such pairs can thus only be destroyed by the 
direct interaction with a third vortex, and not solely by the background noise. 

The results of this section suggest that strictly two-body interactions are not sufficient 
to generate collisions. It is thus natural to study the equivalent three-body problem, which 
is the subject of the next section. 



D. Merging time of a three-body system 

Using Eq. (|4j) and Eq. (|35"D, we can generalize our preceding approach to the case of 
three test vortices submitted to the effective noise induced by far away vortices. For the 
first vortex, the effective equation of motion now reads 

dXi 7/12 „ 3/13 . jf . jf . fAA s 

—7- = -r 2 — - T 3 — + dxiTja + 6y x r) b + v 0x , (44) 
at r 12 r 13 

dyi xu x 13 

-rr = r 2— + T 3 — + oxirjb - dyiTia + v 0y , (45) 
at r 12 r 13 

with similar equations for the two other test vortices. As we are studying the merging time 
of like-sign vortices, we take vortex 1 and 2 to be of circulation +r, whereas vortex 3 is left 
unspecified, with circulation ±r. 

As in the preceding section, these effective equations of motion can be rescaled by ex- 
pressing distances in units of R and times in units of r. In these new units, a collision 
between vortices 1 and 2 occurs when the distance between them is less than the scaled 
dimensionless parameter d c = 2r/R. 
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The sign of the third vortex plays a significant role. When it is the same as that of 
the other two, a phenomenon reminiscent of that which was found for two-body collisions 
occurs: vortices do not collide below a certain radius, at least during numerically observable 
times. 

On the contrary, when this third vortex is of the opposite sign, we observe a smooth 
dependence of the merging time as a function of d c , which is fully compatible with the 
functional form (see fig. 7) 

W ~ (46) 

that is a = 1, and thus £ = 1 using the results of section IV. A. This result is in agree- 
ment with our RG calculation, for which we also found £ ~ 1. Note that for r/R large, 
corresponding to the early stage of the actual dynamics, the apparent value of a is of order 
a ~ 2, which is compatible with an apparent decay exponent in the range £ ~ 0.6 — 0.7 (see 
the inset of fig. 7). 



V. PHYSICAL PICTURE 

The present study strongly suggests that for small surface coverage {r/R <C 1), the 
relevant collision mechanism involves three vortices, one having an opposite circulation from 
the other two. A naive picture would be that of a — r) dipole moving at the typical 
velocity T/r encountering an isolated vortex of any sign. The importance of these fast 



traveling pairs was already suggested in It is likely that the three-body collision 



processes need not to strictly occur in the simple way described above, although this does 
provide an evocative physical picture from which to construct an effective kinetic theory. 
Following the interpretation presented above, the collision rate is 

-77 ~ -na P . x nv dip r, (47) 

where n^p. and t>di P . are the dipole density and typical velocity, and the last term is the 
probability per unit time for a dipole to collide with an isolated vortex. In a simple mean- 
field approach, dipoles of typical size r are formed with density 

n dip ~ n x nr 2 , (48) 

their velocity being of order -u^p. ~ T/r. Finally, we obtain 

n x nVr, (49) 

which corresponds to a merging time 

r merg . ~ {nT x nr 2 Y l ~ (50) 

which in the language of section II. C leads to a = 1, and then to £ = 1 (up to logarithmic 
corrections). 
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This simple interpretation reconciles the paradoxical observation that although vortices 
are hyperdiffusive (see section II.C.l and III.B) the observed decay exponent £ is lower 
than that found for diffusive aggregation |16| , |24] , p5| (see also section 1I.C.2). The dynamics 



is slowed down by the requirement of three close vortices for an actual merging to occur, at 
least for small area coverage. Note that using Eq. ([50]) and energy conservation, we can also 
write 

Tmerg. ~ ~^ ■ (51) 

This relation is particularly well obeyed in the experiment described in [[J, LHS and RHS 
terms being respectively 

, ,t2 j.-0.15±0.Q4 

+ 0.57±0.12 _ i LU ±' 1 + 0.55±0.14 /r \ 

Wg. ~* , and — ~ ^5±01 * " ^ 

Thus, although r merg . does not behave linearly with time in the experiment, which is probably 
due to the fact that the scaling regime was not reached, all the relations found between 
physical quantities (D ~ n" 1 / 2 , v — 1 +£/2 and [i = 3 — £/2 in section II. C, and the above 
relation Eq. ([51])) are fully consistent with our kinetic theory. 

Note that Eq. fl49|) was obtained by Pomeau [15j and by one of us |16| ) but with the use 



of highly questionable physical arguments. In JL3J , the kinetic equation was written in the 
form 

r— ~ — n x nr , (53) 

by arguing that after a time of order r the collision probability is simply the geometrical 
overlapping probability rar 2 , although a cross-section argument is definitely required (see 
section II. C). This argument boils down to the assumption that after a time r, the vortex 
positions can be considered to be randomly generated, and collisions happen if two vortices 
overlap. In [I6| , inspired by the theory of diffusive aggregation p4]j25[] , the kinetic equation 



was written in the form 

but with an incorrect expression for the diffusion coefficient D (found constant in ||I6|). In 
fl6fl , the fact that (v 2 ) is essentially constant was correctly used, but the fluctuation time 
was taken as u~ l instead of r. As we have seen in section IV. C, it happens that just before 
a collision, the effective fluctuation time in fact becomes equal to uj" 1 , which makes the 
agreement of Eq. (p>3) with Eq. (H3) rather incidental. 



VI. CONCLUSION 

In this paper, we have introduced a numerical RG procedure which permits very long- 
time simulations of the vortex Kirchhoff dynamics in a two-dimensional decaying fluid. 
Although we recover a short-time regime compatible with a decay exponent of order £ ^ 0.7 
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when the vortex surface coverage is still large and n/n® > 0.2, we ultimately find a long- 
time asymptotic decay with £ m 1. None of these results can be explained by the simple 
kinetic theory of section II. C based on the occurrence of two-body collisions which predicts 
£ = 4/3. The failure of this "naive" kinetic theory could be explained by our claim that 
strictly two-body collisions are irrelevant for small enough vortex surface coverage (section 
IV. C). For collision processes involving two like-sign vortices and a third opposite-sign 
vortex, we found an average merging time T merg . ~ r/(nr 2 ), fully consistent with a decay 
exponent £ = 1 (section IV. C). A simple kinetic theory based on this collision mechanism 
also leads to £ = 1 and predicts that r merg ~ in agreement with the experiment described 
in ||. Our prediction |21j] that the mean square displacement of surviving vortices goes as 
(x 2 ) ~ t 1+ £/ 2 is in good agreement with our RG simulations and with experiments. Moreover, 
the exponent describing the decay of the flight time distribution is related to £ by ji = 3— £/2, 
in perfect agreement with experiment. 

Our work has so far been limited to the study of the dynamics of a population of vortices 
having the same radii at all times. However, it is important to address the question of the 
possible dependence of the decay exponent £ on the form of the radius distribution and/or 
the initial conditions |2(| , in order to verify its possible universality. Thus, it is a motivating 
challenge to generalize our RG approach to the case of a polydisperse assembly of vortices. 
A crucial point would be to correctly specify the radius of the new vortex reinjected after 
each merging event. This and the study of the effective three-body dynamics of different 
size/circulation vortices should be the subject of a future study |27] . 
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FIGURES 
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FIG. 1. Numerical RG simulations for N = 10,20,40,60,80 (resp. 250,150,30,20,10 samples) 
and the best fit of the N = 20 — 80 curves to the functional form n(t)/no = (1 + t/to)~* (with 
£ = 0.992 and to = 69.8). For clarity, the curves have been offset by an arbitrary factor 2. The time 
t is expressed in units of uj~ 1 and the initial total area coverage is 10%, like for most simulations 
presented in this paper. The inset shows that no/n(t) — 1 is reasonably linear with time (N = 60), 
which is consistent with an exponent £ ~ 1. 



17 




18 




1 00 200 300 

t 



FIG. 3. We compare short time RG simulations (dashed curves, N = 60, 20 samples) with 
previous direct simulations including a polydisperse population of vortices (for which the time unit 
has been scaled to that of the RG simulations). The agreement is good and the apparent decay 
exponent is of order £ w 0.64 (extrapolated to £ ~ 0.72 — 0.75). However, the data display a strong 
curvature in a log-log plot (see inset). 
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FIG. 4. The plot of (v 2 }(t) (RG simulations with N = 80; 10 samples; sampling time 
At = 4w _1 , with no time averaging) displays a slow variation of this quantity fully compati- 



ble with the logarithmic correction obtained in Eq. (26). The thick line is a log-linear fit of the 
scatter plot. 
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FIG. 5. Direct simulations maintaining all radii equal (N = 500, 1000, 2000) are compared to 
RG simulations (N = 40). The saturation to the minimum reachable density u/uq = 2/N is clearly 
seen. However, the direct simulations follow the RG simulations on a longer time domain as N 
increases, up to a time ijy for which some samples have already reached the minimum density. The 
long time apparent exponent for N = 2000 is of order £ « 1.1. 
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FIG. 6. We show a typical long-time relative trajectory within the effective two-body dynamics 
of section IV. B. For clarity, the final time is only i en( j = 200cj _1 , although we have checked that 
the vortices never reach a distance less than ci m i n ~ 0.24i?, up to i en( j ~ lO 5 ^;" 1 . The anisotropy 
of the (square-like) trajectory is due to the use of a Kirchhoff interaction adapted to periodic 
boundary conditions. As the distance d between vortices always remains of order R, these effects 
are visible, although they should disappear when d <^ R (see discussion of section III. A). Still, 
note that when vortices are close to each other, their relative quasi-circular trajectory is not really 
affected by the noise due to other vortices, as explained in section IV. C. 
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FIG. 7. We plot the average merging time r, 



of the dimensionless parameter (R/r) 



[nr 



merjj . 
2\-l 



in units of the mean free time as a function 
as obtained within the effective three-body 



dynamics of section IV.D. We find a good linear behavior for small enough surface coverage. As 
explained in the text, this is consistent with a decay exponent £ = 1. For large surface coverage 
(see inset), a log-log plot leads to an apparent exponent a ~ 1.8, which is compatible with a decay 
exponent £ ~ 0.7, in the language of section IV.A. 
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